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Abstract 



We perform an extrapolation from negative to positive flavour 
\ numbers of full QCD lattice estimates of the rf mass. The extrap- 

■ olations are carried out by keeping p and ir masses at fixed values. 

We find an rf - ir splitting which shows a flavour dependence consis- 
^ | tent with the Witten Veneziano formula based on the U{\) anomaly. 

$_i ' The quantitative splitting is consistent with the estimates made in the 

quenched approximation. 
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1 Introduction 



The quenched approximation in lattice QCD appears to reproduce with re- 
markable accuracy most of the experimental pattern of the hadronic spec- 
trum. Unquenched corrections are small when the comparison with the 
quenched case is performed at fixed values of the lattice spacing and of the 
renormalized quark mass. These two parameters can be fixed by keeping two 
independent physical quantities like the p and 7r masses fixed. The bermion 
method exploits the smoothness of dynamical flavour dependence studied at 
fixed renormalized quantities to make estimates in full QCD which are ex- 
trapolations of results obtained at negative flavour numbers, where fermions 
are replaced by "bermions", i.e. bosons with a fermion action [|l|]. This 
procedure has allowed to show the presence of sizeable unquenching effects 
in light and heavy pseudoscalar decay constants [§]. A quantity where un- 
quenching effects are expected to be dramatic is the rj' mass: it is degenerate 
with the pion mass in the quenched approximation but splits from it in the 
unquenched case. According to the Witten Veneziano formula 

m% =ml + 3jLm* (1) 

where rif is the number of dynamical fermion flavours, N c is the number of 
color and m is related to the topological susceptibility and constant in the 
limit jjl- — > H [§]. For negative flavour numbers the splitting is expected 
to make the rf lighter than the pion. 

In this paper we study the rif dependence of the rf mass with the 
bermion method. We find that eq. [I] is verified: the splitting is indepen- 
dent from rif within the statistical errors and is consistent with the indirect 
estimate made in the quenched approximation, which instead assumes the 
validity of eq. |l]. The correlations of two flavour singlet currents involve two 
classes of diagrams: those where two quark propagators connect the two cur- 
rents (the "fish" diagrams) and those which need gluons to be connected (the 
"eyes" diagrams) . In the quenched approximation only the fish diagrams are 
relevant and they lead to degenerate masses for singlet and non singlet cur- 
rents. The eyes can be evaluated in the approximation where there are no 
intermediate quark loops also in the quenched approximation: this was done 
in [Q by using a modified wall source method. The calculation of the singlet 
correlations with the bermions provides directly eyes and fishes. The first 
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part of this paper is devoted to the detailed discusion of the correlations used 
in the calculations while the second presents our results. 

2 The bermion currents 

The action for QCD with negative flavour number rif can be written in terms 
of rib = 1^/1/2 commuting spinors (f>b(x) called bermions: 

n b 

si nb) [U,<f>,tf] = S G [U] + £ £ri(z)Q(s,z)Q(*,y)&(i/) (2) 

x,y,z i=l 

where Sg[U] is the gauge action, spin and color indices are implicit and Q is 
the Wilson discretization of the Dirac operator 

1 1 3 

[Q4>i]{x] = —l5<Pi{x) - -75 U^(x)(l - 7 M )0i(a: + fi) 

"^75 E Ul(x - a*)(1 + ^)Ux - fi) (3) 

The presence of bermions in the functional integral and in the Monte 
Carlo simulation allows to measure their correlations directly as averages over 
the corresponding fields. In addition, different flavours allow to disentagle 
fish diagrams from eyes and to recombine them in the mixture appropriate 
for singlet currents correlations. 

We have considered the following bermion currents 

Ql\x) = ri(*)7fefc(x) ( 4 ) 

p?(x) = E<lt(v)Q(y>*)M*) (5) 

y 

where % and j are flavour indices. The correlation (P^(x)P^(y)) corresponds 
to the usual fermion correlation with flavour diagonal pseudoscalar currents 
and receives contributions from fish and eye diagrams. Fish diagrams can be 
isolated by measuring the correlation of flavour off-diagonal currents 

G«{x,y) = (Pi j {x)Pt{y)) i*j (6) 
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and eye diagrams by measuring 

H(x,y) = (P*\x)Pi\y)) i / j (7) 
The combination which projects onto the singlet state is: 

G v '(x, y) = Gn(x, y) - n f H(x, y) (8) 

The same correlations with P§ replaced by the pointlike bilinear Q§ of 
eq. f|get the fermion propagators replaced by bermion propagators, which are 
the inverse of the square of the Dirac operator Q 2 . Following the technique 
discussed in || these correlations can be shown to be dominated by the same 
lowest lying states as the P5 ones. In Appendix A we show that the only 
pseudoscalar one particle states which contribute to the correlation 

H Q 2(x,y) = (Q i i(x)Qi j (y)) i^j (9) 

are flavour singlet mesons, and therefore this correlation is dominated at 
large times by the rf without contaminations from pion states. 

In general it is convenient to calculate fish correlations by standard 
inversion of the Dirac operator with fixed origin, while eyes are obtained 
only by direct computation of the correlations between bermion bilinears 
summed over all possible origins. In order to improve the signal of bermion 
correlations, we have applied different variations of the "one link integral" 
technique to the bermion field, replacing, whenever possible, a field with its 
average in the surrounding bermion and gauge configuration. In Appendix 
B we discuss briefly the method and the relevant formulae. 

3 Results 

Our simulations are performed on a 16 3 x 32 lattice at fixed renormalized 
7r and p masses: m n = 0.46 and m p = 0.66 or m w = 0.56 and m p = 0.71. 
We refer to these two cases from the value of the ratio R2 = m^/m 2 which 
is equal to 0.5 and 0.6 respectively: the smaller its value the closer is the 
simulation to the classical limit. 

The main results of this paper are the evidence for a state lighter than 
the pion, shown in figure 1 for R2 = 0.5 and nj = —4, and the determination 
of the parameter of eq. [1| at different values of R2 and as a function of 
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the flavour number. This allows to investigate on further flavour dependence 
of the 77' mass beyond the one already present in the large N c limit in the 
Witten Veneziano expression. 

The determination of the mass splitting Am = m v t — m n can be ex- 
tracted from the quantity 



R(t) = -— log[G v ,(t)/G n (t)] = -— log[l - n f H(t)/G w (t)] (10) 

which grows linearly with time in the region <C t <C ~, where T is the total 
lattice size in the time direction: 

, , Am 

R(t) — > 1 11 

n f 

In the quenched case, by taking the limit rif — > and by assuming 
the validity of Witten Veneziano expression, the previous formulae reduce to 
those used in ||: 

, x Hit) m„ / x 

HM = GW) - 2±J (12) 

Figures 2 to 5 show R(t) as a function of time for R 2 = 0.5 and 0.6 
and different flavour numbers (n/ = and n/ = —4). The straight lines 
are the behaviour expected with the central value of the splitting that we 
extract from our data, fitting R(t) to a straight line for 3 < t < t max with 
5 < tmax < 9. Measurements were taken in a Monte Carlo story consisting 
of the following basic update procedure, 1 heat bath and 3 overrelaxation 
steps for the gauge part and 1 heat bath and 10 overrelaxation steps for 
the bermion part, repeated for a number of sweeps ranging from 23000 to 
141000. The errors on the data points are obtained with a standard jackknife 
algorithm after a binning into clusters of 500 - 1000 measurements. For the 
pure gauge case, in order to keep the advantage of direct field correlations, we 
have introduced pseudofermion fields which are thermalized only in a frozen 
gauge configuration to provide a Monte Carlo inversion of the propagator. 
In order to reach a sufficient accuracy, the auxiliary system is thermalized 
for 200 configurations and the correlation measured for 500 more. It has 
proven useful to introduce many (typically 4) copies of pseudofermion fields 
in order to increase the number of possible flavour recombinations in the 
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correlations of currents P§ and Q l £ . The auxiliary fields are introduced also 
in the one bermion case (n/ = —2) to provide additional non dynamical 
flavours which disentangle eyes from fish diagrams. The results are reported 
in tables 1 and 2 together with some details on the simulation parameters. 
Figure 6 and 7 summarize our results and shows m as a function of the 
flavour number at the two values of criticality that we have explored. The 
results are consistent with an independence from n/ of mo and confirm the 
validity of the expression of eq. 0]. The value of m at R2 = 0.5 is also 
consistent with another determination obtained in the quenched case ||. 
The authors of ref. |J have explored the value of m at values of quark mass 
lower than ours. Their extrapolation to the chiral limit leads to a 7/ mass 
lower than the experimental value by roughly 10%. 

As a further check of the overall picture of singlet meson masses, we 
have also stuidied the singlet and non singlet vector meson correlations. In 
figure 8 the quantity 

R v (t) = -—log[G u {t)/G p (t)] (13) 
n f 

which is the analogous for vector mesons of R(t), is shown for nj = —4 
and i?2 = 0.5: it is compatible with zero, indicating that the two states are 
practically degenerate and confirms the unicity of the pseudoscalar singlet 
state in QCD. 

A last remark concerns the validity of the bermion approach which can 
be monitored through the value of the 7/ mass: at large bermion numbers 
the particle naively speaking becomes massles and may spoil the smoothness 
of the extrapolations from negative flavour numbers. In all calculations that 
we have performed, we have used the rj' mass as a monitor of the limits of 
validity of the method. 
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4 Appendix A 

The bermion action is 

n b 

st b) [U, 0t, 0] = 5 G [C/] + E E (b)t Q 2 (6) (14) 

x k=l 

where and 0^ are commuting fields. For some formal derivations it is useful 
to introduce an equivalent action 

\n f \ 

S ( f - M) [U^,^\ = S G [U]+J2J2^ f) Hxh 5 Q(x,y)^ f \y) (15) 

x,y f=i 

It differs from the standard QCD action only because the ^' and ip fields 
are commuting variables. We consider in the following a theory with at least 
two different bermion flavours (rif < —4). 
The correlation 

H Q2 (t x ) = J d 3 x(Qf(x)Ql>(0)) i ± j (16) 
can be written, by integrating out the bermion fields, as 

H Q2 (t x ) = Z- 1 J VU Tr[ l5 Q-*]Tr[ l5 Q 00 2 ]e- s effM (17) 

where 

Z = j D[U\ e~ s S/V] ( 18 ) 

and 

si n f f[U] = S G [U] + ^Tr(logQ 2 [f/]) (19) 

Using the relation between the bermion propagator Q~y and the fermion 
propagator M~ y l 

^^EtsM-SM^ 1 (20) 

z 

the correlation Hq2(t x ) can be written in the form 

Z- 1 J d 3 x J d'yd'z J VU TriM-^M-^TriM^sM^e-^'M (21) 
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This expression in terms of pseudofermionic fields ipi appears as a four 
point Green function with two of the four variables integrated over the whole 
space-time: 

Jd 3 x Jd 4 yd 4 z(0\T[n 12 (y)a 21 (0)n u (z)a 43 (x)} |0> 
1 Vis = $ipj 



where the ipi fields entering the correlation are labelled with a flavour index. 
The distinction between flavours allows only the contractions represented in 
eq. [T7|. Each new flavour has the same mass, a standard propagator M -1 , 
but does not give rise to additional fermion loops and in general it must be 



considered quenched. However, if \nf\ > 4, the flavours in eq. |22] can be 



associated with the dynamical ip like fields, without introducing additional 
quenched auxiliary flavours. 

The final expression for the continuum case in a infinite volume is ob- 
tained by summing over all different time orderings, by inserting complete 
sets of states in the matrix elements and by using translation invariance. 



For example the contribution to the integrals in eq. 22 from the region with 



< t y < t z < t x can be written as follows: 

EaM I <?x I d 4 yd A ze(t y )e(t z - t y )6(t x - t z ) 
(0|a43(x)|a,p a )(a,p a |7r 3 4(z)|6,p 6 )(6,p b |7ri2(y)|c,p c )(c,p :c |o-2i(0)|0) 

By using translation invariance and by integrating on spatial coordinates this 
expression reduces to 

/ dtydt Z 6(ty)6(t Z - ty)6(t X ~ t Z ) E„ )6 , C ^ ^ % > ~ ^" ^ ^ X 

x<0| < 74 3 (0)|a)(a|7r3 4 (0)|6)(6|7r la (0)|c><c| < 7 21 (0)|0> 

where the sum runs only on states with zero spatial momentum. By inte- 
grating on t y and t z we obtain 

f e ~ mctx —e~ matx e — m b tx —e~ matx \ 

Z^ia,b,c \ (m c -m a )(m c -m b ) (m c -m b )(m b -m a ) J 

x <0|o-43(0) |a> <a|7T34(0) |6> <6|tt 12 (0) |c> (c|a 2 i(0)|0> 

The final result for positive t x , obtained by considering separately twelve 
different time ordered contributions, is 



/ d 3 x J d 4 yd 4 z(0\T [n 12 (y)a 21 (0)7i 34 (z)a 43 (x)} |0) 
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x(0|(743(0)|a)(a|7r 3 4(0)|6)(6|7r 12 (0)| C ) <c|cr 21 (0)|0> 

— (0|7r 34 (0)|a)(a|cr43(0)|6)(6|7ri 2 (0)|c)(c|ff2i(0)|0) 

m a (m c - m b ) 



m c (rrib — m 



(0|a 4 3(0)|a)(a|7r3 4 (0)|6)(6|a 21 (0)|c)(c|7r 12 (0)|0) 



m a m b 
e 



m c m b 

„-m b t 



- <0|7r 34 (0)|a)(a|7r 12 (0)|6)<6| C T 4 3(0)|c){c| O - 21 (0)|0> 
rib 

— <0|flr43(0)|a)(a| < r 2 i(0)|6)(6|7r 3 4(0)|c)(c|7ri 2 (0)|0) 
(0 1 7T 3 4 (0) | a> <a J cr 43 (0) 1 6} <6 1 CT 21 (0) 1 C> <C 1 7T 12 (0) 1 0> 



m a m c 

+ {7T 34 7Ti 2 } (23) 

The main difference with respect to the usual correlations is that the parity 
of the intermediate states is not determined by the one of the operator chosen 
for the correlation: both parity states contribute because of the insertion of 
the parity commuter operator ip^^ip. For large times the lowest lying state 
dominates. In order to determine which states dominate the correlation at 
large times we consider a subgroup SU(2) x SU(2) of the flavour symmetry 
SU(\rif\). One 577(2) subgroup is related to isospin trasformations on the 
flavours 1 and 2 while the other SU (2) is the isospin subgroup for flavours 3 
and 4. The states which can give non zero contributions to the matrix element 
products in eq. |23|, classified with two isospin numbers and parity, are (1,0)^ 
(O,!) 1 * 1 , (0,0) ± , (1,1) ± , including both flavoured pions, corresponding to (1,0)~ 
and (0,1)~, and the (0,0)~ flavour singlet pseudoscalar rj'. However a careful 
examination of eq. ^ shows that the only states whose masses can appear 
in the exponentially decaying factors are flavour non singlet scalar mesons 
(1,0) + and (0,1) + , flavour singlet pseudoscalar (0,0)~ and finally pseudoscalar 
states (1,1) ~ with at least two mesons, while flavoured pions never enter in the 
exponentials. The conclusion is that, if the rj' is lighter then scalar mesons, 
as is always the case in our simulations, the correlation Hqi is dominated at 
large times by this flavour singlet state without contamination from nearest 
states: the pions. 
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5 Appendix B 

We apply to the bermions an extension of the "one link integral" technique 
[0] to derive new bermionic bilinears, which have the same average value but 
smaller variance with the respect to the original ones: the meson correlations 
expressed in terms of these improved bilinears, are found to be correspond- 
ingly much more stable. 

We define the following two points Green functions: 



G Qa {x,y) = (Q$(x)Qg(y)) i?j (24) 
G n (x,y) = (P^(x)Pi l (y)) i / ., (25) 

built from the bilinears defined in eqs. f| and [5[ By integrating over the 
bermionic fields one obtains: 

G Q2 (x,y) = Z~ l !VU Tr[Q-\y,x) l5 Q~ 2 (x,y) l5 }e- s effM (26) 
G n (x,y) = Z- l jVU Tr[Q- 1 {y,x)Q-\x,y)}e- s ^ (27) 



The correlation G n is equal to the standard pion correlation. In the study of 
flavour singlet pseudoscalar states (the rf meson) it is also useful to consider 
the correlations H(x,y) and Hq2(x,y) defined in eqs |7| and ^. H(x,y) can 
be written, after integrating on the bermionic variables, as 

H{x, y) = Z- 1 J VU Tr[Q-\x, x)]Tr[Q-\y, y)] e - s ^^ (28) 

and it isolates the eye diagram. 

We now list the general substitution rules to obtain the bermionic equiv- 
alent of a fermionic bilinear operator. The first is 

i)^(f) and ft^^Qj 5 (29) 
which gives for example the operator P${x). The alternative substitution rule 

4> -> Q<j) and ft -> ft^ 5 (30) 
gives an equivalent operator 

pl?{x) = Y,<t>\(*)Q(^y)Uy) (3i) 
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Given the bilinears 

P 5 (x)=fe^* (32) 
Qs{x) = 4>llh<Px (33) 
we define their improved partners^ where T = 75: 



P^x) = a- 2 [ft( l5 Q - aI)] x T[(rf S Q - al)^} x + — (34) 

a a 

Q' 5 (x) = A- 2 (Q^ - A$)tj(Q 2 $ - A$) x + ^ E r - (35) 



P^'(x) = a'^( l5 Q - aI) 2 } x T[( l5 Q - al) 2 i>] x + ^ £ T aa (36) 



gg(a;) = -a-\q 2 § - A$)tr(g 4 $ - 2AQ 2 $ + a 2 $) x + A- 3 Tr[g^r] (37) 

where 

a -i (38 » 

and 

„ 1 + 16/t 2 , . 

^ = -5?- < 39 > 

The improved operators can replace the original ones with some restric- 
tions on the relative time distance: 

W*)> = TO*)> = ( p 5 (*)> (40) 

(^(^(y)) = (P 5 (x)P 5 (y)) for |t„ - t x \ > 2 (41) 

(Pf(*)if(</)> = (P 5 (x)P 5 (y)) for |t, - t x \ > 4 (42) 

and 

(Q±= ((&(?)) = (Q'&)) (43) 

1 These definitions are valid for a general bilinear with a spinor matrix T. 
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(Q' 5 (x)Q' 5 (y)) = (Q 5 (x)Q 5 (y)) for \t y -t x \>2 (44) 

(Ql{x)Ql(y)) = (Q 5 (x)Q 5 (y)) for \t y - t x \ > 3 (45) 

The above relations are obtained by suitable integration of ip or fields 
in single lattice point which leads to their replacement with a function of the 
surrounding fields. 
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n f 


Tlconf 


/? 


K 


m w 


m p 


R2 


-6 


23 


6.785 


0.156 


0.477(8) 


0.67(1) 


0.507(10) 


-4 


47 


6.463 


0.158 


0.464(4) 


0.658(8) 


0.497(9) 


-2 


89 
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0.161 


0.467(2) 


0.662(5) 


0.498(5) 





92 


5.7 


0.165 


0.457(3) 


0.659(8) 


0.481(7) 


-8 


30 


6.99 


0.153 


0.563(5) 


0.712(10) 


0.625(8) 


-6 


28 


6.71 


0.1548 


0.559(6) 


0.716(9) 


0.610(8) 


-4 


45 


6.4 


0.157 


0.552(3) 


0.710(6) 


0.604(6) 





31 


5.7 


0.163 


0.562(2) 


0.719(4) 


0.611(9) 



Table 1: The parameters of the simulations, the matched values of it and p 
masses and the number of gauge configurations n con f on which meson corre- 
lations are calculated by Dirac operator inversion. 



n f 




P 


K 


m v ,-mir 
n f 


m 


-6 


23000 


6.785 


0.156 


0.035(7) 


0.163(12) 


-4 


121000 


6.463 


0.158 


0.036(5) 


0.168(8) 


-2 


85x500 


6.1 


0.161 


0.029(8) 


0.160(22) 





90x500 


5.7 


0.165 


0.031(7) 


0.167(19) 


-8 


72000 


6.99 


0.153 


0.0146(31) 


0.121(12) 


-6 


28000 


6.71 


0.1548 


0.0155(60) 


0.126(23) 


-4 


141000 


6.4 


0.157 


0.0167(38) 


0.131(14) 





94x500 


5.7 


0.163 


0.0137(36) 


0.124(19) 



Table 2: Results for the rj' — n splitting extracted from the ratio R(t). 
The total number of measurements of the eye correlation H(t) is n meas . For 
rif = and nf = — 2 this correlation was measured 500 times on each fixed 

gauge configuration. In the quenched case we define lim nf ^ mv ' m7T = 
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Figure 1: x: Pion effective mass from the correlation G n for n/ = —4 and 
R 2 = 0.5 (x); flavour singlet pseudoscalar effective mass from the correlation 
Hq2 in the same simulation (*). 
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Figure 2: The ratio R(t) for n f = -4; R 2 = 0.5. 



16 



0.4 



"i i i i i r 



0.35 



0.3 



0.25 



R {t) 0.2 



0.15 



0.1 



0.05 



: 



J I I I I L 



6 7 



1 2 3 4 5 

t 

Figure 3: The ratio R (t) for n f = 0; R 2 = 0.5. 
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Figure 4: The ratio R(t) for n f = -4; R 2 = 0.6. 
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Figure 5: The ratio R (t) for n f = 0; R 2 = 0.6. 
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Figure 6: m as a function of rif for R 2 = 0.5. 
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Figure 7: m as a function of n/ for R 2 = 0.6. 
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Figure 8: The ratio R v (t) for n f = -4; R 2 = 0.5. 
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